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The tight relationship between the masses of black holes and galarxy spheroids in nearby galaxies^ 
implies a causal connection between the growth of these two components. Optically luminous quasars 
host the most prodigious accreting black holes in the Universe and can account for ,1^ 30 % of the 
total cosmological black-hole growth^’^. As typical quasars are not, however, undergoing intense star 
formation and already host massive black holes (> 10® there must have been an earlier pre¬ 

quasar phase when these black holes grew (mass range « 10^—10® Mq). The likely signature of this 
earlier stage is simultaneous black-hole growth and star formation in distant (i.e., 2 > 1; > 8 billion light 
years away) luminous galaxies. Here we report ultra-deep X-ray observations of distant star-forming 
galaxies that are bright at submillimetre wavelengths. We find that the black holes in these galaxies 
are growing almost continuously throughout periods of intense star formation. This activity appears 
to be more tightly associated with these galarxies than any other coeval galarxy populations. We show 
that the black-hole growth from these galaixies is consistent with that expected for the pre-quasar 
phase. 


The most intense sites of star formation at high 
redshift are associated with submillimetre galaxies 
(SMGs)®“®. SMGs are amongst the most bolometrically 
luminous galaxies in the Universe (Lbol ~ 10^® Lq) and 
the majority of the population is coeval with the peak in 
quasar activity (i.e., z 1.5-3.0)”^’®. The apparent associ¬ 

ation of SMGs with some quasars and the similarity in the 
comoving space densities of SMGs and optically luminous 
(Mb < —24) quasars (when corrected for probable source 
lifetimes) provides direct evidence for an evolutionary con¬ 
nection between SMGs and quasars®’^®’®®. However, in 
contrast to quasars, the bolometric output of SMGs ap¬ 
pears to be dominated by powerful star-formation ac¬ 
tivity and any black-hole accretion [i.e., Active Galac¬ 
tic Nuclei (AGN) activity] is comparatively weak^^’^®. 
Given the large available molecular gas supply (typi¬ 
cally Ri 3 X 10^° M©)^^, SMGs can fuel this luminous 
star-formation activity for « 10® yrs (ref 14,15). Vari¬ 
ous pieces of complementary observational support have 
also shown that SMGs are massive galaxies (typically 
« 10^® Mq), suggesting that they will become M, 
spheroid-dominated galaxies in the local Universe®^d6,i7 
AGN activity has been identified in many SMGs®’®^, 
However, the difficulty in obtaining high-quality optical 
spectra and reliable source redshifts has hindered a com¬ 
plete census of AGN activity in SMGs. We have initiated 


a project investigating the properties of AGNs in SMGs 
using deep optical spectroscopic data (obtained with the 
10-m Keck telescope)® and ultra-deep X-ray observations 
(the 2 Ms Ghandra Deep Field-North; CDF-N)®®. As 
X-ray emission appears to be an ubiquitous property of 
AGNs that is comparatively impervious to obscuration, 
the combination of these deep datasets provides a detailed 
census of AGN activity in SMGs. Here we focus on the 
role of black-hole accretion in SMGs and the growth of 
black holes in massive galaxies. The assumed cosmology 
is Mo = 65 km s“® Mpc“®, Dm = 1/3, and Ha = 2/3. 

The parent SMG sample comprises 20 submillimetre 
detected (with the SGUBA camera on the James Clerk 
Maxwell Telescope) sources in the CDF-N region whose 
positions have been reliably located from their radio emis¬ 
sion and which could then be spectroscopically targeted®. 
These sources have overall properties consistent with 
the general SMG population (S'gso/im ~ 4-12 mJy and 
z = 0.56-2.91, with the majority a,t z > 2)®“®. Seventeen 
(« 85%) of the 20 SMGs are detected at X-ray energies. A 
detailed investigation of the nature of the X-ray emission 
has revealed AGN activity in at least 15 (« 75%) SMGs®®; 
the properties of these sources are given in Table 1. The 
X-ray properties of the other 5 (« 25%) SMGs are con¬ 
sistent with those expected from luminous star-formation 
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Figure 1. Black-hole mass accretion rates. SMGs (filled 
circles), and a comparison sample of Mb < — 24 quasars 
(dots) are shown^. The error bars represent the estimated 
uncertainty in the X-ray to bolometric correction for AGN 
activity in SMGs; see Table 1. The solid and dashed 
curves indicate the median and interquartile ranges for 
the comparison quasar sample. The dotted line indicates 
the approximate Eddington-limited mass accretion rate 
for an « 10® Mq black hole (MBH(Edd))- The SMGs have 
mass accretion rates approximately an order of magnitude 
lower than those of coeval quasars^, suggesting that their 
black holes are smaller (typically ~ lO’^-lO® Mq); see text 
for additional black-hole mass constraints. 

activity, but we note that at least one shows evidence for 
AGN activity at near-infrared wavelengths^®. 

Two different observing modes were used in construct¬ 
ing the SMG sample: 14 of the radio sources were specif¬ 
ically targeted with SGUBA observations and 6 of the 
radio sources were detected in blank-field SCUBA maps. 
Although targeting known radio sources with SCUBA ob¬ 
servations could potentially bias our AGN fraction owing 
to contributions to the radio emission from AGN activ¬ 
ity, a statistical analysis does not reveal a strong AGN 
bias in our sample (P=1.0, two-sided Fisher’s exact test); 
see Table 1 for the observing modes of the X-ray classi¬ 
fied AGNs. However, due to the large positional uncer¬ 
tainties of SCUBA sources (w 6-7 arcseconds), redshifts 
could not be obtained for the « 35-50% of the radio- 
undetected SMG population, some of which might host 
AGN activity®’®. Making the conservative assumption 
that none of the radio-undetected SMGs hosts AGN ac¬ 
tivity, our X-ray data alone suggest an AGN fraction in 
the SMG population of > 38tio%- We note that a similar 


fraction of comparably luminous galaxies in the local Uni¬ 
verse host AGN activity^^, although their space density is 
approximately 3 orders of magnitude smaller than SMGs. 

The large AGN fraction in the SMG population indi¬ 
cates that their black holes are growing almost continu¬ 
ously throughout the intense star-formation phase of these 
galaxies [i.e., assuming that all SMGs are growing their 
black holes, the AGN duty cycle is > 381 ^ 0 %]. This sug¬ 
gests that the black holes and galaxy spheroids are grow¬ 
ing concordantly in SMGs^. Such a close association be¬ 
tween AGN and star-formation activity is not seen in the 
coeval field galaxy population or other coeval star-forming 
galaxy populations (w 3-15% AGN fraction)^®’^^. The al¬ 
most continuous black-hole growth in SMGs suggests that 
there is an abundance of available fuel, hinting that the 
accretion may be occurring efficiently (i.e., at or close to 
the Eddington limit), as predicted by theoretical studies 
of the growth of massive black holes®®“®^. The majority 
(« 85%) of the AGNs are obscured (see Table 1), as also 
predicted for black holes undergoing efficient growth^®. 

The unobscured X-ray luminosities of the AGNs (Lx ~ 
erg s“®) suggest that the mass accretion rates 
are modest ( 1 Mq yr“®); see Table 1 and Figure 1. By 

comparison to coeval quasars, the mass accretion rates in 
SMGs are approximately an order of magnitude lower, 
also suggesting modest black-hole masses (i.e., « 10 ^- 
10® Mq; see Table 1 and Figure 1). Although these black- 
hole mass estimates are uncertain, the narrowness of the 
broad emission lines (BELs; typical full-width half maxi¬ 
mum velocities of « 2500 km s“®)®®, when they are visible, 
also suggest black-hole masses of < 10® Mq (ref 4). As 
well-studied AGNs with narrow BELs are found to be ac¬ 
creting at (or even beyond) the Eddington limit'^, this pro¬ 
vides further support for the idea that the mass accretion 
rate in SMGs is approximately Eddington limited. Under 
the assumption of Eddington-limited accretion, the black 
holes in the SMGs will grow by up to an order of mag¬ 
nitude over their ss 10® yr star-formation lifetime. These 
overall properties are consistent with those expected for 
massive galaxies undergoing rapid black-hole growth. 

The total black-hole growth from SMGs can be cal¬ 
culated by integrating the average accretion density over 
the SMG redshift range; see Figure 2. Over the redshift 
interval z = 1.8-3.0 (corresponding to 80% of the SMGs 
studied here), the black-hole density produced by SMGs 
is w 6 I 4 ® X 10® Mq Mpc“®. To put this quantity into 
context we need to compare it to the black-hole growth 
from coeval quasar activity. Considering only quasars in 
the luminosity range of Mb = -24 to -27 (the majority of 
which reside in ^ M* spheroid-dominated galaxies^® and 
have the properties consistent with being the progeny of 
SMGs®), the total SMG black-hole density is « 13^9^% of 
the black-hole growth from quasar activity®. These con¬ 
straints should be considered lower limits since further 
pre-quasar growth at z = 1.8-3.0 may be produced by 
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radio-undetected SMGs and galaxies with fainter submil¬ 
limetre fluxes, if they host AGN activity. However, these 
results are consistent with, for example, SMGs undergo¬ 
ing an intense black-hole growth phase, where the black 
hole is grown from « 10^ to w 10® M©, prior to a high 
accretion-rate quasar phase, where the black hole is grown 
from w 10® to « 8 X 10® Mq. In this scenario, although 
SMGs do not produce a large fraction of the cosmologi¬ 
cal black-hole density, they are responsible for the crucial 
pre-quasar phase where the black holes in massive galax¬ 
ies are rapidly grown. Indeed, the total black-hole growth 
from SMGs could only exceed that from quasar activity if 
the quasar lifetime is insufficient to double the mass of the 
black hole (i.e., < 3 x 10^ yrs assuming Eddington-limited 
accretion). Overall, this picture is in good agreement with 
direct theoretical predictions of the black-hole growth of 
SMGs and quasars^®’^"^. 

What was the catalyst for the rapid black hole and 
stellar growth seen in SMGs? Rest-frame ultra-violet im¬ 
ages taken by the Hubble Space Telescope have shown 
that a considerably larger fraction of SMGs are in ma¬ 
jor mergers (i.e., the merging of two galaxies of compa¬ 
rable masses) than has been found in the coeval galaxy 
population^^’^®. Hydrodynamical simulations have shown 
that major mergers can efficiently transport material to¬ 
ward the central regions of galaxies, providing an effective 
mechanism to trigger nuclear star formation and fuel the 
black hole®°. The result of these major mergers is thought 
to be massive spheroid-dominated galaxies. Ultra-deep X- 
ray observations of SMGs undergoing major mergers have 
shown that AGN activity can be ongoing in both galactic 
components'^. Presumably in these major-merger events 
the black holes will eventually coalesce, further increasing 
the mass of the black hole in the resultant galaxy. 
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Figure 2. Gosmological black-hole accretion density. 
SMGs (filled circle), and Mb = -24 to -27 quasars (solid 
line)^ are shown. Error bars in the x-axis direction in¬ 
dicate the redshift range and error bars in the y-axis di¬ 
rection represent Icr Poisson counting uncertainties. The 
shaded region indicates the uncertainty in the mass ac¬ 
cretion rate conversion for the SMGs; see Table 1. The 
mass accretion density from SMGs is ~ 13lg^% of that 
produced by coeval Mb = -24 to -27 quasar activity over 
the redshift interval z = 1.8-3.0. 
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Table 1. The Properties of the X-ray Classified AGNs 


SMG name 

'S'sSOpm 

(mJy) 

Z 

log(Tx) 
(erg s"i) 

log(LBOL) 

(Le) 

dM 

dt 

{Mq yr-i) 

.^BH(Edd) 

(10" Mq) 

X-ray 

Obscuration? 

Observation 

Mode? 

123549.4 -h621536 

8.3±2.5 

2.20 

44.0 

12.9 

0.281°;?^ 

1 0+2.6 
-‘-•'^-0.9 

Y 

T 

123555.1 -b620901 

5.4±1.9 

1.88 

44.4 

13.2 

0.80lJ;®° 

3.611? 

Y 

T 

123606.7 -b621550 

4.4±1.4 

2.42 

43.8 

12.6 

0.211^1^ 

1.0+^ 

N 

M 

123606.8 -h621021 

11.6±3.5 

2.51 

43.7 

13.1 

A 1 c+0.31 
'^•-‘-^-0.10 

‘ - 0.5 

Y 

T 

123616.1 -b621513 

5.8±1.1 

2.58 

43.7 

13.0 

'^•-‘-'^-0.11 

0 7+^.® 

‘ - 0.5 

Y 

T 

123622.6 -b621629 

7.7±1.3 

2.47 

44.0 

13.1 

A oc+0.70 
U.OO_o.23 

1.6l?;? 

Y 

T 

123629.1 -b621045 

5.0±1.3 

1.01 

43.2 

12.1 

U.UO_o.o3 


Y 

T 

123632.6 -b620800 

5.5±1.3 

1.99 

43.9 

12.9 

0 24+®-^’^ 
'^•^^-0.16 

1 1 + 2.2 

Y 

M 

123635.5 -h621424 

5.5±1.4 

2.01 

44.0 

12.9 

0.301°;®° 

1.4+^ 

Y 

T 

123636.7 -b621156 

7.0±2.1 

0.56 

42.7 

11.1 

0.021^0? 

o.i+l 

N 

M 

123707.2 -h621408 

4.7±1.5 

2.48 

43.8 

12.9 

0.181°;?° 

0.8+? 

Y 

T 

123711.9 -b621325 

4.2±1.4 

1.99 

43.6 

12.7 

n 1 Q+0.25 
'^•-‘-’^-0.08 

U.O-Q 4 

Y 

T 

123712.0 -h621212 

8.0±1.8 

2.91 

43.4 

12.7 


0 4 + 0.8 

0'. + -0.3 

Y 

M 

123716.0 -b620323 

5.3±1.7 

2.04 

44.3 

13.0 

0 ^7+1-33 
'^•'3 ^-0.44 

O.U_2,o 

Y 

T 

123721.8 -b621035 

12.0±3.9 

0.98 

43.7 

11.7 

'^•-‘-'^-0.10 

O'. ‘ - 0.5 

Y 

M 


Notes. — SMG name, submillimetre flux density (SssOMm), redshift (z), and total bolometric luminosity (Lbol) are taken from ref 
9. Unobscured rest-frame 0.5—8.0 keV luminosity (Lx), and the presence of X-ray obscuration (JVh > 10^^ cm“^) are taken from ref 
19. The observation mode refers to whether the source was specifically targeted with a SCUBA observation (“T”) or whether the source 
was detected in a blank-field SCUBA map (“M”): 10 of the ACNs were identified in targeted SCUBA observations and 5 of the ACNs 
were identified in blank-field SCUBA maps. The mass accretion rate (dM/dt) is estimated by converting the X-ray luminosity to the 
ACN bolometric luminosity, under the assumption that the X-ray luminosity accounts for 6^4^% of the bolometric luminosity of the 
ACN (i.e., the range in conversion factors found in ref 20). Although our assumed bolometric conversion factors were originally derived 
from unobscured ACNs, while the majority of the ACNs in our sample are obscured, the important issue is the relationship between the 
intrinsic X-ray emission and the ultraviolet-optical accretion disk emission, which is unlikely to be obscuration dependent. We also note 
that our bolometric conversion factors agree with those estimated for obscured ACNs with similar luminosities (and therefore similar mass 
accretion rates) to the sources studied here®’^^. When converting from the bolometric luminosity to the accreted mass we assumed the 
canonical efficiency of 10% (e = 0.1). The black-hole mass is calculated from the mass accretion rate under the assumption 

that the accretion is Eddington limited; see text for justification. The black-hole masses will be higher if the accretion is sub Eddington; 
see text for additional black-hole mass constraints. 
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